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Comparison of secondary islands in collisional reconnection to Hall reconnection 
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Large-scale resistive Hall-magnetohydrodynamic (Hall-MHD) simulations of the transition from 
Sweet-Parker (collisional) to Hall (collisionless) magnetic reconnection are presented, the first to 
separate effects of secondary islands from collisionless effects. Three main results are described. 
There exists a regime in which secondary islands occur without collisionless effects when the thickness 
of the dissipation regions exceed ion gyroscales. The reconnection rate with secondary islands is 
faster than Sweet-Parker but significantly slower than Hall reconnection. This implies that secondary 
islands are not the cause of the fastest reconnection rates. Because Hall reconnection is much faster, 
its onset causes the ejection of secondary islands from the vicinity of the X-line. These results imply 
that most of the energy release occurs during Hall reconnection. Coronal applications are discussed. 



Magnetic reconnection is widely regarded as the mech- 
anism underlying energy release in the solar corona dur- 
ing flares, coronal mass ejections, and coronal jets 
The Sweet-Parker model [HQ was the first self-consistent 
theory, but is far too slow to explain observations. Much 
has been invested in faster reconnection scenarios, such 
as collisionless (Hall) reconnection 0] in which the Hall 
term plays a key role [H, Hall reconnection seems 
fast enough to explain observed energy release rates 
Q- Lately, the role of secondary islands (plasmoids) on 
Sweet-Parker reconnection has generated much interest. 
While they were discussed some time ago , system- 

atic studies were not carried out until recently. It has 
been argued in various contexts that secondary islands 
make reconnection faster than Sweet-Parker reconnection 
d, IIlTfl3j . (Note, we are discussing secondary islands 
occurring during collisional reconnection, not those that 
occur after collisionless reconnection has begun [HJ|.) 

Understanding secondary islands in Sweet-Parker re- 
connection is important for explaining coronal evolution. 
On the theoretical side, the reconnection rate places con- 
straints on the dynamics. For example, if secondary 
islands make Sweet-Parker reconnection much faster or 
hasten the transition to fast reconnection, it cannot take 
place during pre-flare energy storage. If it remains slow, 
then it can occur while energy accumulates [IM3- On 
the observational side, it was hypothesized that high den- 
sity blobs in current sheets during solar eruptions are sec- 
ondary islands [H, [Hj]. Also, numerous observations of 
reconnection processes display a slow phase preceding an 
eruptive event with an abrupt transition. Examples in- 
clude non-eruptive flux emergence [20| and flows during 
coronal implosion as a result of an impulsive flare pll l22j . 



Past work on secondary islands showed that they ap- 
pear spontaneously due to a secondary tearing instability 
when the Lundquist number S — AncALsp/rjc 2 exceeds 
~ 10 4 10], where L$p is the half-length of the Sweet- 
Parker dissipation region, r\ is the resistivity, and ca is the 
Alfven speed based on the reconnecting magnetic field. 
Equivalently, this can be written as 8 / L$p < 0.01, where 
8 is the thickness of the dissipation region. A study of 
the linear phase of the instability [23[ found a growth 
rate faster than the Alfven transit time along the sheet. 



Recent simulations addressed the nonlinear reconnection 
rate E for high S, showing that it is considerably faster 
than the Sweet-Parker rate and its dependence on S be- 
comes weak [24|, [25j]. However, the simulations only go 
up to S ~ 10 6 , so E is only one order of magnitude faster 
than the Sweet-Parker rate and it is not clear whether it 
will be fast or slow at larger S. Other relevant studies 
showed that E increases with the square root of the num- 
ber of islands [25j, 126J and that secondary islands are sup- 
pressed when reconnection is embedded, meaning that 
the upstream field is smaller than the asymptotic field 
[27| . Many studies consider secondary islands caused by 
external random magnetic perturbations [28l - [3lj |. Other 
studies include the interaction of multiple islands 13 31 and 
a statistical model of multiple island interaction [33[ . 

In addition to increasing the reconnection rate, sec- 
ondary islands hasten the transition to Hall reconnection 
[l^Hj]. When a secondary island forms, the fragmented 
current sheet is shorter, so its Sweet-Parker thickness is 
smaller [2^, [26| . When the layer reaches io n gy roscales 
35, 36], Hall reconnection begins abruptly [H, H 13- 
|39(. This was recently verified using collisional particle- 
in-cell (PIC) simulations HE Ell- 

The only previous study to include both secondary is- 
lands and the Hall effect utilized large PIC simulations 
[IH E3 , but numerical constraints forced S to be small 
enough that Hall reconnection began as soon as a sec- 
ondary island formed. Since the Hall effect arises only 
at ion gyroscales, there should be a regime in which sec- 
ondary islands are present without the Hall effect playing 
a role if the sheet is thicker than ion gyroscales. The goal 
of this study is to separate the two effects and ascertain 
which one leads to dramatically larger reconnection rates, 
which dictates the mechanism releasing the majority of 
the energy during the eruptive phase of reconnection. 

In this paper, the first simulations to separate the two 
effects are presented. There are three main results: (1) 
there is a regime in which secondary islands occur with- 
out collisionless effects entering, (2) the reconnection rate 
is faster than Sweet-Parker but significantly slower than 
Hall reconnection, which shows that secondary islands 
are not the cause of the highest reconnection rates, and 
(3) the onset of Hall reconnection ejects secondary is- 
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lands in the vicinity of the X-line. The latter two imply 
that the majority of the energy release occurs at Hall 
reconnection sites. Coronal applications are discussed. 

Numerical simulations are performed using the two- 
fluid code F3D Magnetic fields and densities are 
normalized to arbitrary values Bq and no, velocities to 
the Alfven speed cao = B$/ (Airrriino) 1 / 2 where rrii is the 
ion mass, lengths to the ion inertial length dio — c/co P i, 
times to the ion cyclotron time fi~- , electric fields to 
Eq = caqBq/c, and resistivities to 770 = 4ircAodio/c 2 . 

The initial configuration is a double tearing mode 
with two Harris sheets, B x (y) = tanh[(y + L y /4)/wo} — 
tanh[(y — L y /A)/wo\ — 1, where wq is the initial current 
layer thickness and L y is the system size in the inflow di- 
rection. Total pressure is balanced initially using a non- 
uniform density which asymptotes to 1. The tempera- 
ture T = 1 is constant and uniform. A single X-line is 
seeded using a coherent magnetic perturbation of ampli- 
tude 1.6 x 10 -2 to rapidly achieve nonlinear reconnection. 
Initial random magnetic perturbations break symmetry 
so secondary islands are ejected. There is no initial out- 
of-plane (guide) magnetic field. Boundaries in both direc- 
tions are periodic. Electron inertia is m e — mi/25. This 
value is acceptable since we focus on the onset of Hall 
reconnection at ion scales rather than electron scales. 

Simulation parameters are chosen so reconnection will 
proceed in three distinct phases: Sweet-Parker with- 
out secondary islands, Sweet-Parker with secondary is- 
lands, and Hall reconnection. A very large system size 
L x x L y = 819.2 x 409.6 is employed with resistivity 
r) = 0.008, corresponding to a global Lundquist number 
S g = L x /r) ~ 10 5 , which exceeds the Biskamp criterion 
of 10 4 . To postpone secondary island onset, we choose 
Wq = 12.0 which makes the reconnection embedded [27| . 
Embedding makes the Sweet-Parker layer thicker since 
S ~ {rjLsp / CAup) 1 / 2 , where CAup is the Alfven speed 
based on the upstream magnetic field B up . For wide cur- 
rent layers, B up ~ Bq8/wq [27| . so eliminating B up gives 
8 ~ (t]L sp w ) 1/3 ~ 2.7, where L S p ~ L x /4 ~ 200 m our 
periodic system. Thus, the layer begins wider than d%, 
and since S/Lsp > 0.01, no secondary islands occur ini- 
tially and the system undergoes Sweet-Parker reconnec- 
tion. The reconnection inflow convects in stronger mag- 
netic fields, so the current sheet self-consistently thins. 
Islands arise when S/Lsp ~ 0.01, which gives 8 ~ 2.0 
since L S p ~ 200. It has been shown [H, S that if N 
X-lines are present, 8 decreases by a factor of TV 1 / 2 . For 
a single secondary island (N — 2), the layer shrinks to 
5 ~ 2.0/2 1 / 2 ~ 1.4. This exceeds d u so Sweet-Parker 
with secondary islands should persist. Hall reconnection 
only starts when 8 ~ 1, so three distinct phases occur. 

A simulation is first performed with a grid scale A = 
0.2 and the results are qualitatively consistent with ex- 
pectations. To assure A does not play a role, the simu- 
lations are redone with A = 0.1, giving comparable re- 
sults. Data is presented only from the high resolution 
runs. The equations employ fourth order diffusion with 
coefficient D4 = 1.75 x 10 ~ 4 to damp noise at the grid 
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FIG. 1: (Color) (a) Reconnection rate E as a function of time 
t. The solid (blue) line is a Hall-MHD run. Dashed lines 
at t ~ 700 and 1780 indicate the onset of secondary islands 
and Hall reconnection, respectively. The dot-dashed (red) line 
shows E for a simulation restarted at t — 1120 with no Hall 
effect and m e = 0. (b) Thickness 8 of the dissipation region 
vs. t. Horizontal dotted lines mark predicted 8 for the onset 
of secondary islands (8 ~ 2) and Hall reconnection (8 ~ 1). 



scale. A smaller value of D4 leads to a slightly larger Hall 
reconnection rate, but does not alter our key conclusions. 

We now summarize the simulation results, followed by 
a careful justification of the conclusions. At early times, 
Sweet-Parker reconnection prevails. A secondary island 
first appears at t ~ 700. Reconnection proceeds with 
the secondary island until t ~ 1780, when Hall reconnec- 
tion onsets. Thus, reconnection proceeds in three dis- 
tinct phases including an extended phase with secondary 
islands but without the Hall effect triggered. 

We compare the reconnection rate E in the three 
phases to each other and to theoretical predictions in 
Fig. HJa), showing E vs. time t as the solid (blue) line. 
We measure E as the time rate of change in the difference 
in magnetic flux function ip between the main X-line and 
O-line. Dashed lines at t = 700 and 1780 denote where 
secondary islands and the Hall effect arise, respectively. 
Ignoring secondary islands, Sweet-Parker theory predicts 
E - E SP - 0.006, where E SP - {^/L SP ) 1/2 , and L SP ~ 
200. This assumes the magnetic field is its asymptotic 
value of 1. The measured value is E ~ 0.004, slightly 
lower than predicted as expected because B up < 1 (the 
reconnection is embedded). When TV X-lines are present, 
E scales as E S i ~ E S pVN [H, The measured rate 
of 0.005 is consistent with this for a single secondary is- 
land (N = 2). After the Hall effect onsets, E increases 
by about an order of magnitude. Therefore, the recon- 
nection rate with secondary islands is faster than Sweet- 
Parker, but significantly slower than Hall reconnection. 

The transitions occur when predicted, as shown in 



3 




FIG. 2: (Color) Time history plot of the out-of-plane current FIG. 3: (Color) Time evolution of J z from the other current 
density J z in the outflow direction. Dashed lines mark when sheet in our double tearing mode setup, showing the ejection 
a secondary island appears and when the Hall term onsets. of secondary islands when Hall reconnection onsets. 



Fig.QJb). We plot 8, measured as the half- width at half- 
max of J z in the inflow direction through the X-line, vs. t. 
The dotted lines at S ~ 2 and 1 show the predicted value 
when islands and the Hall effect should appear, respec- 
tively. These conditions are met at t ~ 700 and 1780, in 
good agreement with the observed transitions. 

The appearance of new physics can be seen in direct 
observations of the out-of-plane current density J z . A 
two-dimensional time history plot of J z in the outflow 
direction is plotted in Fig. [5] Only the half domain cen- 
tered on the seeded X-line is shown. The raw data is 
sampled at a rate of one frame per 70 time units, so lin- 
ear interpolation is used to smooth data between time 
slices. The effect is cosmetic, not substantive. The color 
bar is stretched to enhance visibility of weaker currents. 

Early in time, J z is structureless and extends the half- 
length of the domain, as expected during Sweet-Parker 
reconnection. A secondary island near x = appears 
as a dark spot with associated strengthening of the frag- 
mented current sheets. This occurs at t ~ 700, marked 
by the vertical dashed line. This agrees with Biskamp's 
criterion shown in Fig. [If b). As time evolves, the island 
grows and 5 shrinks. When S ~ di, Hall reconnection 
onsets and the current sheet becomes much shorter and 
intense, appearing as a sharp peak in J z in Fig. [2j This 
begins at t ~ 1780, as also marked in Fig. Hfb). 

There are two locations where Hall reconnection on- 
sets. An X-line near x ~ —70 onsets slightly earlier than 
an X-line at x ~ 70. As Fig. [5] vividly shows, the lat- 
ter X-line is ejected from the dissipation region, along 
with the secondary island, which is ejected at the Alfven 
speed. The ejection of the secondary island implies that 
the two effect will not (locally) coexist, so most of the 
energy is released at Hall reconnection sites. 

This current sheet has only a single secondary island 



and one may ask whether this result remains valid in 
more realistic settings with multiple islands. To address 
this, we show results from the other current sheet in our 
double tearing mode setup, which self-consistently devel- 
ops multiple islands. Figure [3] shows J z at three times 
near the onset of Hall reconnection. Panel (a) is just as 
Hall reconnection onsets at x ~ 20, showing three pre- 
existing secondary islands. The Hall reconnection X-line 
grows steadily, as shown in panel (b). Panel (c) shows 
that the single X-line at x ~ 20 is the only one to persist 
as all of the secondary islands are ejected. This suggests 
that the ejection of nearby secondary islands by Hall re- 
connection sites is a robust result, and may reasonably 
represent local behavior in a macroscopic current sheet. 

A careful determination of when the Hall effect begins 
to become important is obtained using a time history 
plot of the out-of-plane Hall electric field Eh z = J y B x /n 
in the inflow direction through the main X-line, plotted 
in Fig. Ufa). (Note, this cut is in the inflow direction, 
while Fig. [2] is in the outflow direction.) The color bar 
is again stretched. The plot clearly shows that Eh z does 
not contribute during the secondary island phase. A cut 
of Eh z in time, taken at the solid (gray) line in Fig. U{a) , 
is plotted in Fig. HJb). The onset time, defined as when 
Eh z reaches 1% of its maximum value, is at t ~ 1780, 
the time that E begins to increase as seen in Fig. [lja). 

To emphasize differences between Sweet-Parker with 
secondary islands and Hall reconnection, we restart the 
simulation at t — 1120 with the Hall effect and electron 
inertia disabled. The reconnection rate is plotted as the 
dot-dashed (red) line in Fig. HJa). The value reaches 
E ~ 0.009 as the asymptotic upstream field reaches the 
dissipation region, in excellent agreement with the pre- 
dicted value E S i ~ E SP \fN ~ 0.009 with N = 2 for 
a single island. This rate is consistent with the largest 
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FIG. 4: (Color) (a) Time history plot of the out-of-plane Hall 
electric field Ehz in the inflow direction, (b) Plot of Ehz vs. t 
at the y location marked in panel (a). 

scaling studies done to date [24j. Note, E remains nearly 
an order of magnitude slower than Hall reconnection. Al- 
though the present evidence is based on simulations only 
up to S ~ 10 5 , it is clear that secondary island reconnec- 
tion does not produce the fastest reconnection rates. 
The present results may be relevant for observations of 



two-phase reconnection events in the corona. In observa- 
tions of flux emergence [2(| , a slow phase of reconnection 
preceded an abrupt transition to a fast phase ~ 30 times 
faster (compare the slopes in their Fig. 18). In observa- 
tions of the contraction of magnetic loops in an impulsive 
flare the contraction velocity abruptly increased by 
a factor of ~ 16. It is enticing to attribute these obser- 
vations to a transition from resistive secondary island re- 
connection at a normalized reconnection rate of E ~ 0.01 
(consistent with implications of Refs. 0, H^]) to Hall re- 
connection ~ 10 times faster which occurs abruptly when 
gyroscales are reached. The existing level of accuracy of 
both theory and observations make such an identification 
premature, but it remains an exciting possibility. 

Assumptions in this work that require further study 
include using a Spitzer resistivity instead of the uniform 
resistivity employed here, including Ohmic heating and 
viscosity, and including Dreicer field effects, which may 
be important for the transition to Hall reconnection. The 
simulations have no guide field, but one would be ex- 
pected in the corona. Also, the simulations are two- 
dimensional; three-dimensional effects are not included. 
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